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Abstract. Relativistic simulations in 3+1 dimensions typicaUy monitor the 
Hamihonian and momentum constraints during evolution, with significant violations 
of these constraints indicating the presence of instabilities. In this paper we rewrite the 
momentum constraints as first-order evolution equations, and show that the popular 
^ ' BSSN formulation of the Einstein equations explicitly uses the momentum constraints 

O '< as evolution equations. We conjecture that this feature is a key reason for the relative 

$-H . success of the BSSN formulation in numerical relativity. 
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ff^ . 1. Introduction 
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I Three decades after the earhest attempts to study black hole spacetimes numerically, 

I the problem of simulating the collision and coalescence of two black holes remains at the 

P I forefront of numerical relativity. Only in the last year have the first stable, long-term 

numerical simulations of black hole collision and coalescence in three-dimensions been 

X ; achieved [H [2]. The primary focus of numerical relativity has at last shifted to the 
$_( ' 

' astrophysical implications of these events, with recent simulations investigating recoil in 
the final product of binary black hole mergers, and the extraction of gravitational wave 
profiles for comparison with observations [3, HI [5]. 

Despite these advances, progress is still needed in understanding both the 
mathematical and computational formulations of the equations. The most popular 
formulation of the Einstein equations in numerical relativity is the BSSN system [6l [7] 
(but see [8] for an alternate approach). The Einstein equations are split into distinct 
evolution and constraint equations, with the constraints being used to construct initial 
data, after which they are monitored (rather than enforced) as the evolution proceeds. 
In addition, a variety of technical and computational techniques are employed to achieve 
stable evolutions, including causal differencing, artificial diffusion, grid excision, moving 
punctures, special shift conditions, and many others which have been tried over the last 
few decades. 

In this paper we consider the role played by the momentum constraints in evolution, 
and in particular, their appearance in the BSSN formulation. It is well known that the 
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momentum constraints are vital to the stability of the BSSN system, where they are 
used to eliminate the spatial divergence of the trace-free extrinsic curvature from one of 
the evolution equations. This is a necessary condition for stability [7]. The Hamiltonian 
and momentum constraints are then monitored during evolution, and give an indication 
of the quality and stability of the simulation. 

The close link between unstable simulations and uncontrolled violation of the 
constraints has led to the development of several constrained evolution algorithms 
[HI [ini [m Il2]- These schemes are designed to enforce the Hamiltonian and/or 
momentum constraints throughout the evolution, either through direct projection or 
gentle numerical relaxation onto the constraint surface. 

We will show that the BSSN system is already a partially constrained formulation of 
the Einstein equations, since the evolution equations used for the conformal connection 
functions are the momentum constraints rewritten as first-order evolution equations. 
It is in this sense that the BSSN formulation is a partially constrained evolution 
scheme, with a form of the momentum constraints being actively enforced throughout 
the evolution. 

We begin by briefly reviewing the BSSN formulation and its derivation, highlighting 
the role played by the momentum constraints. In section 3 the momentum constraints 
are recast as first-order evolution equations for a particular combination of the metric 
and its spatial derivatives. Then, in section 4, a conformal transformation on this 
variable demonstrates that the momentum constraints may be used to evolve the 
conformal connection functions. A more direct calculation demonstrates that the 
resulting equation is identical to to the standard BSSN evolution equation for T\ In 
the final section we consider the consequences of this observation. 

2. The ADM and BSSN evolution equations 

The Baumgarte-Shapiro-Shibata-Nakamura (BSSN) formulation of the Einstein 
equations [6l [7] can be derived from the canonical ADM form of the equations. With 
the metric in the 3-1-1 form 

ds"^ = -(a^ _ p.pi^ + 2[3i dt dx' + gtj dx'dx\ (1) 

the components of the three-metric gij become the dynamical variables, and the 
remaining coordinate freedoms are expressed by the lapse function a and shift vector 
The Einstein equations consist of six second order (in both time and space), non- 
linear partial differential equations, together with an additional four constraints which 
only contain first temporal derivatives of the gij. By introducing the extrinsic curvature 
tensor Kij the Einstein equations can be written as a set of twelve evolution equations 
which are explicitly first-order in time (neglecting source terms for simplicity), 

{dt - Cp) gij = - 2aKij (2) 
(dt - Cp) K,, = a [Rij - 2KiiK'^ + KK^^) - WiV.a, (3) 
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where Cp is the Lie derivative along /3* and K = g^^Kij. The constraint equations are 

R+{Ky -K,jK'' = 2p, (4) 
VjK"^ - g'^VjK = S\ (5) 

and are known as the Hamiltonian and momentum constraints, respectively. These 
are mathematically conserved by the evolution equations, in the sense that if they are 
satisfied on some spacelike hypersurface then solutions of the evolution equations will 
continue to satisfy the constraints for all time. 

This 3+1 form of the Einstein equations is known as the ADM system, and was 
the basis of initial efforts to construct numerical solutions. Until quite recently all 
attempts to evolve black holes in three dimensions were unsuccessful (see [I3] for a 
review) , with the impending failure of the code signposted by rapid growth in violations 
of the constraints. This led many investigators to propose and test new formulations 
of the equations, among the most popular being the so-called BSSN system developed 
by Shibata and Nakamura [6], and later popularized by Baumgarte and Shapiro [7]. 
Their approach splits the extrinsic curvature tensor into its trace (K) and trace-free 
(Aij) parts, and then performs a conformal decomposition of the three-geometry, with 
9ij — ^^^9ij ^iid Aij = e^'^ Aij. Then 

Ki, = e^^ (A, + l~g,,K), (6) 
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and for consistency we require that det(^jj) = 1 and A\ = g^^Aij = 0. 

The ADM evolution equations ([2]) and (I3l) are transformed to evolve gij and Aij, 

{dt - Cp) (jij = - 2aAij (7) 

{dt - Cfs) A, = e-^^ [aR,, - WiVjd]^^ + a[K Aj - 2 A^i^ (8) 

where "TF" denotes the trace-free portion of the bracketed expression. Similarly, the 
ADM equations lead to 

{dt-Cp)(t) = - ^-aK (9) 

{dt -Cp)K = a [a,jA^ + ^ - (10) 

which are used to evolve (j) and K . 

The BSSN approach also introduces the conformal connection functions 

r = -d,~g'\ (11) 

for which evolution equations can be obtained by differentiating with respect to the time 
coordinate and commuting the partial derivatives. This yields 

dtV = -dj {2aA' + C0^) (12) 

where equation ([7]) has been used to eliminate explicit time derivatives of the conformal 
three-metric. The Lie derivative is given by 

= P'dk~g'' - 2f^'dk(3^^ + ^ -g'^dkP', (13) 
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since g^^ is a tensor density of weight |. 

In principle this completes the system of evolution equations, since the conformal 
connection coefficients could be evolved using equation ( |T2l) . However, it is found 
that the system is only more stable than the original ADM formulation if the spatial 
divergence of A^^ is eliminated from equation (JT2l) . This is achieved using the momentum 
constraints [7]. Under the conformal transformation described above, the momentum 
constraints take the form 

djA'^ + 6A'^dj(P - \ -g'^djK + V.^A^^ = e^'*'S\ (14) 

where 5"* is a source term. Using these to eliminate the spatial divergence of A^^ from 
(IT2D gives 

dtV = - 2A'Wja + 2a (gA'^ d^cj) - ^ ~g'^ d^K + f}^ A^^ - S'^ (15) 

- d, (^P'dk~g'' - 2f^'dkP'^ + ^ -g'^dkP' 

The standard BSSN evolution equations consist of (|7l)- (nil together with (HM . and 
are used to evolve the BSSN variables gij, Aij, (j), K and F*. This system has been found 
to be relatively stable, and is basis for most of the successful recent binary black hole 
calculations [Hill [5]. 

3. The momentum constraints as first-order evolution equations 

The momentum constraints ([5]) are a set of three differential constraints relating the 
three-metric and extrinsic curvature. Since the extrinsic curvature is itself related 
to the time development of the three-metric through its definition, equation ([2]), the 
momentum constraints can be viewed as equations involving the time derivative of gij. 
In this section we explicitly rewrite the momentum constraints as first-order evolution 
equations for a function of the three-metric and its spatial derivatives. 

The momentum constraints with source terms, equation ([5]), can be written as 

djK'^ + V]^K''' + Vi^^K'^ - g'^djK = S\ (16) 

where T^jj^ are the Christoffel symbols and represents the sources. Using the definition 
of extrinsic curvature, equation ([2]), the first term in the constraints becomes 

djK^' = d, (dtg^^ - Lf^g'^ 
and by commuting the spatial and temporal partial derivatives of gij we find that 

Applying a similar procedure to the gradient of K, we have 

g'^d.K = ^ (-a, (g'^g^'^gab) - 2Kg'^d,a + g'^'K^'d^g,, 
+ g'^g'^'d, {Cpgab) + g'^ d.gabCpg'^' + g'^'d.gaM^ 
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Combining these results we can write the momentum constraints fll6l) as 

dt {d,g'' + (f 'd, \ng) =2 (k^^ - g'^ k) d,a 

+ 2a [K'^dj \ng + S' - Vi^K'^ - r}^K^'=) 

+ d, {Cpg'^) + g'^d, {g^'Cpgab) + (d, \ng) L^g^^ 

where we have used g'^^djgah = dj In g and g = det{gij). 

The momentum constraints can thus be viewed as the natural evolution equations 
for the subsidiary variable 

A' = d,g'^+g'W,\ng, 

the first term of which is suggestive of the conformal connection coefficients defined 
by equation (ITT]) . The momentum constraints written in their most general form as 
evolution equations are then 

dtA' = 2 (k'^ - g'^K) d.ja + 2a [k'^Oj In g + S'- Vi^K'^ - T)^K^^) 

+ d, [Cpg^^) + g^^d, {g^^'L^gah) + {d, \ng) ll^g^\ (17) 

where A* is a subsidiary variable since it is not truly independent, but rather a function 
of the metric and its spatial derivatives. We conjecture that using the momentum 
constraints to evolve A* in the ADM formulation, alongside the evolution of gij and Kjj, 
would enforce the consistency of the constrained evolution system more thoroughly than 
evolving the three-metric and extrinsic curvature alone. 

Finally, we note that equation (ITTIl is essentially the same form of the momentum 
constraint used in the first-order flux-conservative formulation of Bona and Masso |14j . 
The Bona-Masso approach introduces a new variable which is evolved using the 
momentum constraint, and in our notation Vi = \ gijAK 

4. Momentum constraints in the BSSN formulation 

In the previous section we demonstrated that the momentum constraints can be used to 
evolve the subsidiary variable A*, which is itself a function of the metric and its spatial 
derivatives. In this section we show that under a conformal transformation A* is closely 
related to the BSSN variable F*. Guided by this insight, we will derive the evolution 
equation for F* directly from the conformally decomposed momentum constraints. 

Beginning with the variable A* introduced above, the conformal decomposition used 
in section 2 gives 

which suggests that we define A* = e'^'^A*, so that 
A« = _f ^ - 8 ~g'^ dj(l). 

It is clear that A* splits into the BSSN conformal connection coefficient and the gradient 
of the conformal factor. 
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It follows from the results of the previous section that the momentum constraint 
can be used to evolve F*. By calculating the time derivative of A* we find that 

and eliminating the dtcj) and dtg^^ terms, using equations (JTj) and ([9]), relates the time 
evolution of A* directly to the evolution of F*. Substituting this into the momentum 
constraints in the form of equation (ITTI) will then give an evolution equation for F*. 

This shows, somewhat indirectly, that the momentum constraints can be used to 
evolve the conformal connection coefficients F*. It is more enlightening, however, if we 
begin with the conformal form of the momentum constraints, 

djA'^ + QA'^dj<i) - ^ -g'^djK + f)^A^'' = e^'^S' (18) 

and perform a calculation analogous to that of the previous section. 
It follows directly from the definition of Aij, equation ([7]), that 
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This allows us to re-express the first term in the momentum constraints as 

d^A^' = ^( - d,{Cpr) ) - ^ A^'d^a, 

and commuting the partial derivatives of g'^^ and using the definition of F* we can write 

djA'= = ( dr + 2i'^a,a + d,{Cpr) ) . 

Using this result to replace the divergence of A^^ in the conformal momentum constraints, 
equation (ITSil . we find 



- ^{dtV + 2A'^d,a + d,{Cp~g'' 

+ QA'^dj(t) - ^ ~g'^djK + f]j^A^^ = e^^S\ 

It is important to note that these are precisely the momentum constraints; all that we 
have done is replace A^^ with its definition in terms of g'^K Finally, expanding the Lie 
derivative and rearranging yields the momentum constraints in a new form: 

dtV = - 2A'^dja + 2a [qA'^ d^cj) - ^ ~g'^ diK + f }^ A^^ - S'^ 

- dj (^P''dk~g'' - 2f^%(3^^ + ^ -g^^dkP" 

This is identical to equation f|T5l) . the standard equation used to evolve the 
conformal connection coefficients in the BSSN formulation. From this we conclude that 
the usual statement that the BSSN formulation is stable when the spatial divergence of 
A^^ is eliminated from the evolution equation for F' [71 is equivalent to the statement 
that the BSSN system is stable when the momentum constraints themselves are used to 
evolve the conformal connection coefficients. 
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5. Conclusions 

We have shown that the BSSN formulation utilizes the momentum constraints as 
evolution equations for the conformal connection functions F*. It follows that the 
momentum constraints should be satisfied to computational accuracy throughout the 
evolution, although in practice violations of the "definitional constraints" 

r + dj~g'^ = 0, det(^,,) = l, and g'^A,=0, (19) 

could lead to errors when the momentum constraints (|T8l) are calculated during 
evolution. 

We conjecture that the BSSN formulation shows superior stability properties to the 
ADM system because errors in the definitional constraints do not correspond directly 
to spurious sources of momentum being pumped into the simulation. Although the 
Hamiltonian and momentum constraints are mathematically conserved by the ADM 
evolution equations, small computational errors prevent the strict consistency of the 
system. The particular advantage of the BSSN system is that the conformal connection 
functions provide a way of more strictly and actively enforcing the consistency of the 
equations, provided the violation of the Hamiltonian constraint is not too great. 

We expect that in practice any violation observed when calculating the momentum 
constraints (flSj) during evolution arises from inconsistencies introduce by violations of 
the new definitional constraints ( |T9|) . It is possible that additional stability in the system 
may be gained by including the addition terms in the evolution equation for which 
"drop out" , but may in practise be non-zero. For example, terms containing 

g^^djcjab = djlng 

are nominally zero because g = det(^jj) = 1, but in practise we expect small violations 
in one of more of the constraints ( fT9|) . 

Our result is consistent with recent work by Marronetti [HI [12] , where a constraint 
relaxation method is applied to the BSSN formulation. Although weak enforcement of 
the Hamiltonian constraint was found to provide significant improvement in the results 
[TT] . there was little change when a similar approach was applied to the momentum 
constraints [12]. As we have seen, the momentum constraints are already being 
directly enforced in the BSSN approach. It is likely that the small improvements seen 
in Marronetti's momentum relaxation scheme were due to improved (albeit indirect) 
enforcement of the remaining constraints F* + djg^^ = and det(^jj) = 1. 

We have shown that the BSSN formulation ensures that the momentum constraints 
are enforced throughout the evolution. The generalized harmonic approach of Pretorius 
[8] is also, by construction, constraint damping. Since all of the successful codes 
currently used in numerical relativity are based on one of these formulations, it appears 
that active constraint enforcement is a vital ingredient for the success of any new 
formulation of the Einstein equations designed for use in numerical relativity. 
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